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(f) ■ Abstract. - We have considered the interaction of a pair of spatially separated two-level atoms 

with the electromagnetic field in its vacuum state and we have analyzed the amount of entangle- 
ment induced between the two atoms by the non local field fluctuations. This has allowed us to 
characterize the quantum nature of the non local correlations of the electromagnetic field vacuum 
state as well as to link the induced quantum entanglement with Casimir - Polder potentials. 

Ok 
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The zero - point fluctuations of the vacuum state of the electromagnetic field arc char- 
acterized by strong non local correlations [1-3] which are at the origin of phenomena like 
Casimir - Polder forces [4,5]. An interesting open problem is the possibility to characterize 
, the quantum nature of such non local correlations. So far most of the efforts done towards 

this direction have focused on violations of suitable Bell's inequalities by the vacuum state 
fluctuations [6]. However a direct experimental detection of such inequalities violation for 
the vacuum state is awkward. In this letter we shall take a somewhat different approach. 
It is a well known fact that when two quantum subsystems, e.g. two atoms, interact with a 
common bath, they become entangled (see for instance [7]). Such a pair of subsystems can 
therefore be used as a probe of the non local vacuum field fluctuations. In other words the 
quantum nature of such fluctuations can be characterized by the amount of entanglement 
induced between the two spatially separated probe atoms. Some interesting contributions in 
this direction have already appeared in literature [8] . In the following we shall quantify the 
entanglement induced between the two probe atoms by means of the concurrence [9] as this 
is amenable to a straightforward physical interpretation [10]. Indeed we will show that the 
concurrence turns out to be linked to the Casimir - Polder potentials. Casimir-Polder forces 
are long-range interactions between neutral atoms or molecules arising from their interac- 
tion with the common electromagnetic radiation field in its vacuum state. For atoms in the 
ground state the Casimir-Polder potential behaves as i?~ 6 for interatomic distances smaller 
than a characteristic distance of the order of an appropriate average of the atomic transition 
wavelengths - but large enough to neglect any overlap between the electron wavefunctions - 
and as R~ 7 for larger distances. 

Casimir-Polder potentials have been experimentally detected in several physical systems. 
Here we mention two experiments, which are closer to our analysis [11, 12], in which the 
deflection of neutral atoms due to the interaction with conducting plates, has been measured. 
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It is possible to show that the Casimir-Polder potential is also strictly related to the spatial 
correlations of the vacuum fluctuations of the electromagnetic field. In fact, it can be 
obtained from the classical interaction of the instantaneous dipole moments of the two atoms, 
which are induced by the spatially correlated zero-point quantum field fluctuations [13]. 
Furthermore the use of (three-body) dynamical Casimir-Polder forces for investigating the 
nonlocality of field correlations has been also recently suggested [14] 

To introduce all the mathematical tools used in our analysis let us briefly recall the def- 
inition and basic properties of the concurrence. Quantum Information Theory has provided 
new powerful mathematical tools to quantify the amount of entanglement between quantum 
subsystems in some specific situations. In particular for the case of a mixed state of two 
two- level systems the entanglement of formation £ F is a suitable entanglement quantifier [9] , 
as it quantifies the amount of non local resources needed to create a given state. From a 
mathematical viewpoint the main advantage of considering the entanglement of formation 
is the fact that it is a monotone function of the so - called concurrence, a quantity which 
is relatively easy to calculate. For an arbitrary bipartite system described by the density 
operator p, the entanglement of formation £p turns out to be equal to [9] 

£f(p) = -xlog 2 x- (1 - z)log 2 (l - x) (1) 



where x = (1 + \Jl — C 2 (p))/2 and the concurrence C(p) is defined as 

C(p) = max {0, ct\ — a 2 — a% — a^} . (2) 

where {aii} (i = 1, ..,4) are the square roots of the eigenvalues (in non-increasing order) of 
the non-Hermitian operator p = p(a y <8> <J y )p*{<J y (g> a y ), a y is the y-Pauli operator and p* is 
the complex conjugate of p, in the basis of a z operator. Since the entanglement of formation 
is a monotonic function of the concurrence we will use the latter as entanglement quantifier. 
Furthermore, as anticipated, the concurrence turns out to be directly linked to the Casimir 
- Polder potential, and therefore amenable of an experimental detection. 

The system we have considered consists of a pair of spatially separated two-level atoms, 
A and B, placed at and R# respectively, at a distance R = R^ — Rs large enough to 
neglect any overlap between the electron wavefunctions, interacting with the the electromag- 
netic field in its vacuum state. The atom-radiation system is described by the multipolar 
Hamiltonian in the dipole approximation 



with 



H = Hab + Hp + Hab.f (3) 
H AB = huj Q ^ S« (4) 

i=A,B 

H F = y^huka^akj (5) 

kj 

tt ST X^\( c ®n ^k-Ri , ,(*)* t „-ik-Ri c( 

HAB.F - 2^ 2^ [{ e kj a ^3 e S + +%' «kj e S - 



=A,B kj 

\ e -ik.R i5 W +e W* akje ik.R i5 W^ (6) 



where loq = ckg represents the separation in angular frequency of the two atomic levels, S z , 
S + and S- represent the atomic pseudo-spin operators, a^j and denote the annihilation 
and creation operators, respectively, of photons with wavevector k and polarization j, and 
the coupling constant, in the multipolar representation, 
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is purely imaginary for real atomic dipoles and linear polarization unit vectors ekj (note 
that no rotating wave approximation has been made). The two atomic states considered 
might be, for example, the Is ground state and one of the 2p excited states of a hydrogen 
atom, which satisfy the electric dipole selection rules for the emission of a photon. In this 
case, all the three degenerate p excited states equally participate to the interaction with the 
field, of course. However, in the following we shall consider the interaction of the ground 
state with only one of these states, because this makes easier the connection of our results 
for the concurrence with the Casimir-Polder forces. On the other hand, a sum over the three 
degenerate excited states can be done at the end of the calculation, restoring the spherical 
symmetry of the atoms; in fact, it is known from the theory of Casimir-Polder forces that 
the relevant virtual transitions to the excited states occur independently, up to the required 
order in perturbation theory. This makes the two-level system considered here a convenient 
system to obtain and analyze the physically relevant quantities. 

The normalized, dressed ground state of the two atoms can be written at second order 
of approximation in the electric dipole 0{d 2 ) in the compact form 

W) = c ag\ 0S{°kj}) + Y ° e a^ \ e lk J> + Y c 3 e ^ gel ^ + Cee ' ee {° k j}> 

kj kj 

+ \ Y Y ''-w.k.kv I gglkjlk'j') + ^ Y Y c ee,kjk'j' I eelkjlk'j') (8) 

kj k'j' kj k'j' 

where \g) and \e) denote, respectively, the ground and excited state of each atom. The 
explicit form of the various probability amplitudes will be given later on. We remark the 
different properties of the first- and second-order corrections (i.e., the terms with odd and 
even number of photons, respectively): whereas the former give local effects, such as the 
Lamb shift of the energy levels of each atom, the latter are responsible for the Casimir-Polder 
interaction and are thus essentially of non local character. 

The reduced density operator of the two atoms qab = Tr fi e i,i\ip){'ip\, obtained by tracing 
over the field variables, has nonvanishing entries only in the two diagonals of its matrix 
representation. For this class of density operators the concurrence C can be expressed in 
terms of physical quantities, such as correlation functions and average values [15]: 



C = 2max{0,Ci,C 2 }, (9) 




where gij is defined in terms of expectation value of spin operators as 

9ij = (S^Sf), M i = {Sf + Sf)/2, SSi = (Sf-Sf) (12) 

Such quantities can be evaluated with the help of the wave function eq.(8) therefore dispens- 
ing us from an explicit calculation of the density operator qab- Straightforward algebraic 
manipulations give a simpler form for C\ and C2 for our two-atom system: 

Ci = y /(\ee){gg\)(\gg)(^ -y /(\eg)(eg\)(\g^e^ (13) 
C 2 = y/(\eg){ge\)(\ge)(eg\) y/(\ee)(ee\)(\gg)(gg\) (14) 



p-3 



M.A.Ci 



G. Compagno 1, G. M. Palma 2, R. Passantc 1, F. Persico 1 



Our goal is the evaluation of the concurrence at second order of approximation. Since 
(\gg)(gg\) = 1 + 0(d 2 ), the average value (|ee)(ee|) must be evaluated at fourth order of 
approximation. From eq.(8) we obtain 



(\ ee )(gg\) 
(\eg)(ge\) 

(\eg)(eg\) 
(\ee)(ee\) 



(\99)(ee\)* = c gg c* ee , 

(lff e }( e S|}* = J2 c egMl C *ge,kj' 



kj 



I C ep,kj I 7 

kj 



(\9e)( 9 e\) 



El \ C geMj 
kj 



\ c ee\ + „ E/ l C ee. k J' k 'i'l 



k,j k',j' 

The probability amplitudes required for the evaluation of the concurrence are 



c egMj 



Cge.kj — 



£ kj e 

h(aj + 0J k ) ' 

IB) -ik-R B 
% e 

ft(w + Wfe) ' 

E 



1 x - £ kj £ kj l e 



kj 



(^) (B) „-»k-R A -ik' R B 
£ kj £ k'j' e e 

ft 2 (w + Wfc)(wo + Wfe') 



7i(wo + w/s) 

(A — B) 



With the help of these expressions eq.(14) becomes 



(15) 
(16) 

(17) 
(18) 



(19) 
(20) 
(21) 

(22) 



C 2 



U) (B) ik -R M) (B) ik' R 

E £ kj £ kj e e k'j' e k'j' e 
ft 2 (',.,„ _i_ ,.„ ^2 

kj 



1/2 



,(B) p ik'.R 

"kj "kj ^ £ k'j' £ k'j' e 
^ ft 2 (u, + ^ fe ) 2 ^ ^ 2 ( W o+w fe 2 



(A) (B) ik.R 
£ kj e kj e 



E 

kj 



^)|2 



ft 2 (w +Wfc) 2 ^ ^ 2 (W +Wfe/) 2 



C (B) 12 
E k'j' 



+ Ce 



1/2 



Thus C 2 is negative and the concurrence is C = max(0, 2 Ci), with 



E/ Ce 9> k j E/ C 9 e M'j' 
kj k'j' 



1/2 



(23) 



(24) 



The sums in square brackets formally diverge when the sum runs from k = to k — > 
oo; however, as noted above, the probability amplitudes c eg ^j and c ge ^j describe local 
interactions of each atom with the vacuum field and play no role in the interaction between 
the two atoms. Such terms are responsible only for the appearance of local energy shifts 
and do not contribute to interatomic forces nor to non local correlation of the fluctuations 
of atomic dipoles. Indeed, in the research literature on the interaction of pairs of atoms 
with vacuum fluctuations, all such terms which do not depend on the atomic separation 
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(and therefore diverge) are neglected - see e.g. [1, 16, 17] . Following this approach the 
concurrence is 

C = 2d = 2\c ee \ (25) 

Therefore for any separation R the concurrence, and thus the entanglement, is determined 
by c ee only. This fact is consistent with the observation that the term c ee | ee{0kj}} in the 
dressed ground state Eq.(8) is entirely due to the interaction between the two atoms via the 
common electromagnetic field. 

When we transform the sums over the wave vector k into integrals following the usual 
prescription V^ 1 — ► (27r)~ 3 J d 3 k we obtain 



E 



M) (B) ±ik . n 
h(u>Q + LOk) 



7T ' * 



(26) 



where the differential operator 



D 



( V S mn -\- V m Vr, 



l 
R 



(5 mn R m Rn) gR2 + (S mn ZR m R n )(^ R2 RgR j 



(27) 



has been introduced for ease of notation and f(x) denotes one of the auxiliary functions of 
the integral sine and cosine functions [18]. Thus the concurrence takes the form 



C 



irhojo 



J2 d ™ d nD*J(k R) 



(28) 



Using a terminology typical of the Casimir - Polder context, this expression simplifies 
considerably when we examine its behavior in the near zone, defined by koR -C 1, and in 
the far zone, defined by koR S> 1: 



^ |d^d*-3(d^-R)(d*-R)| 

C (nearzonc) ~ > k ° R « 1 ( 29 ) 

C ~ 8c|d*-d*-2(d^R)(d*-R)| 

C (farzonc) - 2^1 > ^ » 1 (30) 

The above expressions, which are the main result of this letter, are amenable of a straight- 
forward physical interpretation. To illustrate the relations between our expressions for the 
concurrence and the Casimir - Polder potential we first observe that this latter can be written 
in the following form [13]: 

W c -P = 3? E tt° k I £ ki( R B)»^( R A)< I Ok)a A (k)a B {k)Vt m (k,R)] J (31) 

\kj,tm I 



where 



2irhc 



(0 k | E^(R B ) m E^(R A )i | fc ) = (e kj ) m (e kj ) e fce* kR (32) 

is the equal-time spatial correlation function of the electrical field modes in the vacuum 
state and evaluated at the position of the two atoms , 

a ^ - W^t) (33) 
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is the dynamic electric polarizability of the atoms and the quantity 

cosfci? 



V tm (k,K) = -D 



R 



R 



= k 



( ~ ~ \ cos kR ( ~ - \ f sin kR cos kR\ 

[S em R e R m ) - [S em ZReRm) [j^ + -j^- ) 



(34) 



is the classical interaction potential between two dipoles oscillating at frequency ck [19]. 
The form (31) of the Casimir-Polder energy emphasizes the role of the spatial correlations 
of the electromagnetic vacuum state. 

The link between the near - zone concurrence and the near zone Casimir - Polder R~ 6 
potential can be easily established by noting that the latter coincides with the well-known 
van der Waals potential between two neutral atoms. This was first derived by London 
[20] and can be obtained treating by second order perturbation theory the dipole - dipole 
interaction hamiltonian [21] 

d^d 5 -3(d^R)(d 5 -R) 

R? ' [6b > 

In the near - zone the Casimir -Polder potential is therefore essentially of electrostatic nature. 
Note that indeed that the expression (29) for the concurrence in the near zone, is the ratio 
between the interaction energy between two permanent dipoles and the energy separation 
hujo- 

The fact that in the near - zone the interaction between the two atoms is essentially 
of electrostatic nature and can be described by an effective Hamiltonian in which the field 
degrees of freedom are eliminated suggests that in such region the state of the two atoms is 
described by a pure density operator. Indeed, in the near zone, we have p 2 AB = pab up to 
second-order approximation if only the leading terms in the expansion in series of powers of 
koR are kept. This implies that the field induces pure bipartite entanglement between the 
probe atoms i.e. there is no entanglement between atoms and field but only entanglement 
between the atoms which is mediated by the field. 

The "far - zone" R~ 7 behavior of the Casimir-Polder potential, eq.(31), stems from 
retardation effects, and it is a typical manifestation of the quantum nature of the electro- 
magnetic field. In an analogous way, a change of the power law (R~ 4 instead of R~ 3 ) is 
found in the concurrence when we move to the far zone (see eq.(30)). Here entanglement 
can be interpreted as a consequence of vacuum correlations of the electric field E [22], since 
the concurrence can be cast in the form 

C (far zone) = 2^ ^"71 ' k ° R » 1 ( 36 ) 

The analogy with the Casimir-Polder potential becomes more transparent if we note that 
eq. (31 ), in the far zone, can be cast in the following form:" 

^jAjA jB jB 

W C-P(far zone) = E " ^ E ^^ j (n A ) n )V ii (k, R) (37) 

mni£ kj 

Both (36) and (37) show clearly the role of the spatial correlations of the quantum field 
vacuum fluctuations in correlating the atomic dipoles. The only difference is that while the 
concurrence (36) is due to spatial correlations between the fields at positions R^ and 
while the Casimir - Polder potential (37) is due to the correlations between the field modes 
at positions Ha and Rg. Note incidentally that in the far zone the system (atoms plus field) 
manifests an essentially tripartite entanglement, since in this region we have p\ B ^ pab- 

In summary we have evaluated the entanglement between two spatially separated two- 
level atoms interacting with the vacuum state of the radiation field. In particular we have 
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shown that the concurrence can be cast in a form with strong analogies with the Casimir 
- Polder potential. The same experimental setup used to measure such potential can be 
used to infer the amount of entanglement between two spatially separated neutral atoms. 
To estimate the expected order of magnitude of the concurrence one can assume that the 
two atoms are two hydrogen atoms. This gives in the near zone Cn.z. ~ (R/ao)~ 3 where 
ao is the Bohr radius, while in the far zone Cf z « a(R/ao)~ 4 where a = e 2 /(hc) is the 
dimcnsionless fine structure constant. 

The above discussion provides a physically transparent characterization of the entangle- 
ment created by the non local zero-point field fluctuations and suggests a strategy for its 
experimental detection. 

* * * 
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